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Abstract 



The FAC, PMS, and BLM optimization methods are applied to the QED 
corrections to the muon lifetime in the Fermi V-A theory. The FAC and PMS 
scales are close to m e , while the BLM scale nearly concides with the geomet- 
ric average ^/m e m^. The optimized expressions are employed to estimate the 
third order coefficient in the a(m^) expansion and the theoretical error of the 
perturbative series. Using arguments based on effective field theory and a 
simple examination of Feynman diagrams, it is shown that, if contributions 
of 0{amj 1 /My i r) are neglected, the corrections to muon decay in the SM fac- 
torize into the QED correction of the Fermi V-A theory and the electroweak 
amplitude <? 2 /(l — Ar), both of which are strictly scale- independent. We use 
the results to clarify how the QED corrections to muon decay and the Fermi 
constant Gf should be used in the SM, and what is the natural choice of 
scales if running couplings are employed. 
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I. INTRODUCTION 



The one-loop radiative corrections to /x-decay in the Fermi theory were evaluated approx- 
imately four decades ago. These studies included the correction to the electron spectrum 
and muon lifetime |E],p| , as well as the momentum dependence of the asymmetry and the 
integrated asymmetry for polarized muons [lj. Although the Fermi theory is, in general, 
non-renormalizable, the /i-decay case is especial. For vector and axial vector interactions of 
the charge retention order, which are the relevant couplings in the two-component theory of 
the neutrino and in the V-A theory, a theorem assures us that, to first order in G^, but all 
orders in a, the corrections to /z-decay are convergent after mass and charge renormalization 
0. A striking cancellation of mass singularities in the correction to the lifetime and inte- 
grated asymmetries was found, which occurs also for the scalar, pseudoscalar, and tensor 
Fermi interactions, and the /3-decay lifetime |l| . These observations were one of the motiva- 
tions in the derivation of the KLN theorem Years later, in order to clarify a controversy 
that arose in the determination of the Fermi constant, the corrections of O {a 2 In {m^/m e )) 
to r M were obtained 0. 

Very recently, in an important theoretical development, van Ritbergen and Stuart com- 
pleted the evaluation of O (a 2 ) corrections to r M in the local V-A theory, in the limit 
m t/ m fi ~> 0U- Their final answer can be expressed succintly as 

C^) = ^c l+ (^)\, (1 ) 

7T y 7T J 



ci = \ (f - * 2 ) ; c 2 = 6.700 . (2) 

In Eq.(JI|), a{m^)c\l'K is the one-loop result |TJ and a(m M ) is a running coupling defined by 

ct 

aim u ) = 7 rr . (3) 

\ 37r 2n z J m e 

In Refs. @,f|, the contribution of the last term in the denominator of Eq.(JD is separated 
out, as (a 3 /27r 2 ) In (m M /m e ). The two expressions differ by only 0.2 ppm, so that we will 
employ the more succint expression of Eq.(|]). 



In Section 2, we apply the traditional FAC 0, PMS ||10|| , and BLM |TT| optimization 
schemes to the expansion of Eq.([|). We use those methods for estimations of the third order 
coefficient and the theoretical error in Eq.([TJ). In Section 3 we turn our attention to the 
radiative corrections to muon decay in the Standard Model (SM). Using arguments based 
on effective field theory and simple considerations of Feynman diagrams, we show that, if 
terms of 0(am 2 l /My ir ) are neglected, the radiative corrections factor out into [1 + C{m e )] 
(defined in Section 2) and the electroweak amplitude g 2 /(l — Ar), which are separately 
scale- independent. We use the results to discuss the application of the Fermi constant Gf 
in electroweak physics. Section 4 presents our conclusions. 



2 



II. APPLICATION OF OPTIMIZATION METHODS 



As it is well known, a(fi) does not run below the m e scale and a(m e ) is identified with 
the conventional fine structure constant, a -1 = 137.03599959(38)(13) fl2| . Setting then 
a = a{m e ) in Eq.(||) and replacing m e — > fi, we have 



aim 



1 - 



3tt 2tt 2 



ln^ 



(4) 



Inserting this expression into Eq.(|l]), expanding and truncating in second order, we find 



7T 



7T 



2 , rria 
c 2 + -d In 

3 /i 



(5) 



The method of Fastest Apparent Convergence (FAC) || chooses the scale /i FAC in such 
a manner that the NLO term vanishes. Thus, we have In {m^/ ^, FAC ) = — 3c2/(2ci), which 
leads to /x FAC = 0.801 m e . As a(/z) does not run below ji = m e , we interpret this result as 
/i FAC = m e . Eq.(|) becomes 



C(a*fac) = C(m e ) = -ci + - 

7T \7T/ 



.2 m M 
c 2 + -ci In 

3 m P . 



(6) 



an expansion in terms of the fine structure constant a. The (9((a;/7r) 2 /n(m /1 /m e )) term 
concides with the expression found in Ref. 0. Recalling m^/m e = 206.768273(24) [f[^| , we 
see that the two terms between square brackets nearly cancel each other and Eq. (§) becomes 



C(/w) =C(m e ) 



a 
— < 

7T 



"Ci + 



0.26724 



(7) 



The Principle of Minimal Sensitivity (PMS) |10| identifies fi PMS with the stationary point 
of C(/x). Applying n(d/djj) to C(/x), and recalling the renormalization group equation (RGE) 



d . . 2 2 a 3 (/x) 

Gt/i 37T 27T Z 



we obtain 



(8) 



n—c(fi) = -( a ^ 

d\i 3 V 7r 



3 2 m„ 

c 2 + o ci + -Ci In 

8 3 n 



(9) 



Thus, the stationary point is given by 

/^PMS 



In- 



m. 



3 /3 £2 

2 V8 ci 



which leads to yU PM s = 1.40636 m e . We note that Eq.® is consistent with 



an expression that can be obtained by replacing — > /i, /i — > m e in Eq.(|]), and will be 
useful in the following discussion. Inserting /i = /i PMS in Eq.(|), we obtain 

n , x _ a(/i PMS )ci /a(/i PMS )\ 2 3 

L'l^PMsJ — qCi, 

7T y 7T y O 

or 

C(/* PMS ) = a(/iFMs)Cl + ( V 0.67868 . (12) 

7T \ 7T / 



The BLM method [[TIJ chooses /i BLM so as to cancel the terms in c 2 proportional to 
rif, the number of light fermions. In the /i-decay case, there is only one light fermion, 
the electron. From Ref. we learn that the electron loop contributions to c 2 (virtual 
loops and pair creation) is 3.22034. Splitting c 2 = 3.22034 + 3.47966, one chooses fi BLM to 
cancel the first contribution in the second order coefficient of Eq.(|3p. Thus, In (m M //i BLM ) = 
-(3/2)(3.22034/ci), which leads to /i BLM = m^/14.4267 = 14.3323 m e . We note that this 
is very close to ^ym e m^ = 14.38 m e , the geometric average of the two masses. The BLM 
expansion is 

COW) = + I ^HA) 2 3.47966 . (13) 

7T \ 7T / 




The FAC, PMS, and BLM optimization expressions of Eqs.([7|, 0, [13|), can be used to 
estimate the coefficient C3 of (a(m M )/7r) 3 in the expansion of Eq.(Jl|). It is well known that 
such estimations are not particularly useful if contributions of an important class of new 
diagrams open up at the relevant level. For instance, in the radiative corrections to g — 2, 
large contributions due to the light by light scattering open up already in O (a 3 ). In /x- 
decay we encounter a more felicitous situation, as light by light scattering contributes only 
in 0(a A ). Writing the three optimized expansions in the generic form 

C(/i*) = ^c 1+ (^] c*, (14) 

and expressing a(/i*) in terms of a(m^) by replacing m M m M in Eq.([|), we find 

the estimation 

M^h^ + ii + i*)^ (15) 

Using the values of /1* and c\ obtained above, Eq.(|T5]) leads to 

f -19.9 (FAC) 

(c 3 ) est = -20.0 (PMS) (16) 
[ -15.7 (BLM) 

We see that the three estimates are quite close. 
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It is also interesting to compare the numerical results of the optimized expansions among 
themselves and with that of Eq.([|). Using the values of fi PMS and /i BLM found before, and 
evaluating a((j,*) (/i* = Li PMS ,Li BLM ) via Eq.(|i~0l), the expansions of Eqs. ([7|, p^ , pr^ ,p]) lead to 



C(m e ) 

BLM J 

C(m„) 



-3 



-3 



-4.202402 x 10 
-4.202403 x 10 
-4.202348 x 10~ 3 
-4.202147 x 10~ 3 



(FAC), 
(PMS), 
(BLM), 



(a(m M J exp. 



(17) 
(18) 
(19) 
(20) 



We see that the three optimization methods give very close results, with a maximum absolute 
difference of 5.4 x 10~ 8 . As the important mass scales for ii decay in the Fermi theory are m e 
and m M , it is natural to consider m e < fi < as the relevant range. We may therefore take 
the maximum difference among the four evaluations above as an estimate of the theoretical 
error. This is given by the difference C(m M ) — C(fi PMS ) = 2.6 x 10~ 7 , which is equivalent 
to the estimated third order contribution 20.0[a(m M )/7r] 3 in the PMS estimate [Cf.Eq. (|T6|) 1. 
This leads to an error of 1.3 x 10~ 7 in the determination of Gp, which is very close to the 
estimate given in Ref. || from the consideration of the known leading third order logarithm 
in C(m e ). We also recall that the one- loop corrections proportional to powers and logarithms 



of m 2 /m 2 can be obtained by combining Ref. jl3| and Ref. H], and have been reported in 



Refs. Hjl4|1. When the tree-level fx decay phase space is factored out, the leading correction 



of this type is 



a m i 
7r m 



2Alog^ - 9 - 4vr 2 



4.3 x 10' 



(21) 



In very high precision calculations, Eq.(^l|) should be added to Eqs. (|17| - ^0D so that the 
expansions, rounded to four decimals, become 



C(m e 



PMS J 
BLM,/ 



-4.1981 x 10" 
-4.1980 x 10 
-4.1978 x 10 



-3 



-3 



(22) 
(23) 
(24) 

In 2 m M /m e ) has been esti- 
In fact, we note that if 



The error due to the uncalculated terms of 0((a/7r) 2 (mg/m 2 t 
mated to be a few times 10~ 7 ||. This seems quite resonable 
the C((a/7r) 2 (mg/m 2 )) contributions had the same magnitude relative to the leading 0(a 2 ) 
term (a(m^)/ir) 2 C2 in Eq.(|I]), as the (9((a/7r)(m 2 /m 2 )) bear with respect to (a(m^)/7r)c 1 , 
their contribution would be even smaller, a few times 10~ 8 . 



III. FACTORIZATION OF THE RADIATIVE CORRECTIONS TO MUON 

DECAY IN THE SM 



In the on-shell renormalization scheme of the SM, it is customary to write 1/r^ in the 
form 
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1 _ P g 4 [1 + C(m e )] 
-Ar) 2 ' 



32 rr^v 



P = f 



mz 



in: 



, 3 m 2 
1 + M 



5M, 2 



m. 



/(x) = l-8a;-12x 2 lna; + 



192tt 3! 



lx 3 - x 4 , 



(25) 

(26) 
(27) 



where g 2 



e 2 / sin 2 



sin 



M^/M 2 Z , M w 



and M 2 are pole masses, C(m e ) is the 



radiative correction of the Fermi V-A theory and Ar is the electroweak radiative correction 
introduced in Ref. [fLER . Before the work of Refs. only the terms involving c\ in Eq.@ 
were known and C(m e ) was approximated by the expression 



a 

-ci 

7T 



l + 2a ln m, 
37r m P . 



In Eq. (|26|) , /(m 2 /m 2 )m^/1927r 3 is a phase space factor and 3m 2 /5M^ is the tree-level 
contribution of the W-propagator. 

If terms of C»(am 2 /M^) are neglected in the radiative corrections, and Ar is approx- 
imated by Ar^ — e 4 ReIl2^(M§), where Ar^ is the one-loop contribution to Ar, and 



e 4 H^\Mz) is the two-loop contribution to the renormalized vacuum polarization function 



at q 2 = Mf , it was shown in Ref. |T(| that Eq.Q2l)|) contains all the one-loop corrections, as 
well as all two-loop contributions involving mass singularities (i.e. logarithms mm, where m 
is a generic light fermion mass). It also contains all the terms of O [(a ln(M^/m))™]. In fact, 
if so desired, all these logarithms can be absorbed by expressing g 4 /(l — Ar) 2 in terms of 
the running coupling ol(M%) = a/(l — Act), where Act = — Rell^Mf ). This result follows 
from the observation that, in g 4 /(l — Ar) 2 , such logarithms arise from the renormalization 
of «o in terms of a. 

We now show that, if contributions of 0{am 2 L / M^) are neglected in the radiative cor- 
rections, the factorization displayed in Eq . (|25"D , involving the QED corrections of the Fermi 
theory and the electroweak factor g 4 / (1 — Ar) 2 , is valid to all orders in perturbation theory. 
This applies not only to the corrections to the lifetime, but also to those affecting the elec- 
tron spectrum, as well as the momentum dependence of the asymmmetry and the integrated 
asymmetry in the case of polarized muons. We present two arguments, one based on the 
effective field theory approach (|I1J , the other involving a simple discussion of higher order 
Feynman diagrams. 

If contribution of 0{am 2 l / M^) are neglected in the radiative corrections and the tree- 
level term 3m 2 /5M^ in Eq.(p6|) is for the moment disregarded, the effective field theory at 
the muon mass scale is the local V-A four-fermion Lagrangian density 



C 



Gp 

"71 



^(l-y 5 )^ * 7/i (i 



plus QED, plus QCD. Therefore, one can systematically evaluate the corrections to the 
spectrum, lifetime, and asymmetry in muon decay on the basis of this effective Lagrangian. 
This procedure results in the usual expressions involving Gf and the radiative corrections of 
the Fermi V-A theory |l| . As the latter are convergent to all orders in a || , there is no need 
to cancel ultraviolet divergences and, therefore, in their expression, there is no reference 
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to the high mass scale Mz of the underlying theory. As mentioned in the Introduction, 
these corrections are known at the one-loop level in the case of the electron spectrum and 
asymmetry, and have now been evaluated at the two- loop level in the case of l/r M . In 
particular, one finds 



1 _ G 2 F m 5 

t7 ~ 1927T 3 



F"V tl 2 / 2 



f (mi/mi) [1 + C(m e )\ , (28) 

where C{m e ) is the two-loop correction discussed in Sections 1 and 2. This leads to Eq.(p5[) 
provided we identify 

i I 3 m M 

Gf_ g 2 i + iQM^ (9q] 

V2 8M 2 , (1 - Ar) • 1 yj 

Eq.(p9|) has a very simple interpretation: it is the matching relation that expresses the 
coupling constant of the effective, low energy theory, in terms of the coupling constants 
and radiative corrections of the underlying theory. It is very important to note that C(m e ) 
and g 2 /(l — Ar) are separately scale-independent quantities. The //-independence of C(m e ) 
follows from the fact that, when expressed in terms of physical parameters such as a, m e , m^, 
it is convergent to all orders of perturbation theory ||. The /i-independence of g 2 /{l — Ar) 
follows then from the fact that it can be expressed in terms of physical observables via 
Eqs. (28,29). As explained above, in the on-shell scheme of renormalization g 2 is defined 



in terms of e 2 , M w and M z , which are physical quantities and, therefore, /i- independent. 
It follows that Ar is \i- independent, an important property that has been verified at the 
one-loop level [IS|, and through terms of 0(g A M 2 /M$y) at the two-loop level [18]. 



The same conclusion, concerning the factorization of QED and electroweak corrections 
when terms of 0{am 2 l / M^) are neglected, can be reached by a simple analysis of Feyn- 
man diagrams. Consider, for instance, a two-loop diagram involving a photon of virtual 
momentum k attached to external /i and/or e lines and a second loop that includes heavy 
particles. The relevant momenta in the QED correction of the Fermi theory are \k 2 \ < m 2 . 
If we neglect terms of 0(am 2 l /My V ), we can set k = in the heavy-loop integration and 
the two-loops factor out. We also note that when \k 2 \ becomes a few times larger than 
m 2 , we can set p e = 0, where p e is the electron four-momentum. All reference to p e is lost 
and we see that such domain of virtual momenta does not contribute to the corrections 
to the electron spectrum. Thus, the latter involves essentially the domain |fc 2 | < m 2 , for 
which the factorization is valid. A second consideration, based on Feynman diagrams, refers 
to the factorization of QED corrections relative to the dominant electroweak corrections. 
These arise from the renormalization of the bare coupling constant g 1 = Cq/sq. Consider, 
for simplicity, the one-loop photonic diagrams to //-decay in the SM. Each diagram carries a 
factor eg/sp associated with the virtual W interchange. The renormalization of Cq/sq in such 
diagrams leads, in higher orders, to dominant electroweak corrections that clearly factorize 
with respect to the QED corrections. The same arguments can be extended to higher orders. 

The crucial Eq.fl2"9|) can be simplified by introducing a coupling constant 

G ° r 



11 1 _L 3m M 
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so that Eqs. fl28p9D become 



V2 



G 2 m 5 

9 2 1 
8M 2 , (1 - Ar) ' 



, 3 m 2 
1 + M 



5M^ 



[1 + C(m e )] , 



(30) 
(31) 



These are, in fact, the expressions most commonly used in literature. The difference between 
Gfj, and Gf is 0.5 ppm; it is completely negligible at present and it would only be of marginal 
interest if the experimental error in r M is reduced by a factor 10. Nonetheless, this factor is 
frequently included in the theoretical expressions, as it is a tree-level contribution from the 
SM. Eq. (|30|) has the convenient feature that all the kinematical factors of the SM calculation 
are separated out. 

The arguments carried above can be applied mutatis mutandi to other renormalization 
frameworks such as the MS scheme of Ref. ||19|| , in which couplings in g 2 /(l — Ar) are 
identified with running MS parameters evaluated at Mz, masses are still interpreted as 
pole masses, and Ar is modified accordingly. 

The authors of Ref. || discuss a renormalization framework in which the couplings in 
g 2 /(l — Ar) are also identified with running parameters evaluated at Mz- In particular, 
they claim that, in order to ensure the consistency of the renormalization scheme, the 
QED correction C of the local theory, when applied to the SM, must be evaluated via 
Eq.fll]) with ot(m M ) — > a e (M z ) (a e (M z ) is the expression obtained from Eq.(|]) by replacing 
m M — > Mz)- This claim has a curious consequence: before c<i was evaluated, the current 
experimental value of r M led, via Eq.(f28|), to G F = 1.16639(2) x 10~ 5 /GeV 2 . With the 
incorporation of the C2 term in Eq.([l]) one obtains Gf = 1.16637(1) x 10~ 5 /GeV 2 . However, 
according to the new claim, the correction C is to be evaluated differently when applied to 
the SM than in the Fermi theory, with the result that the relevant value is turned back to 
Gf = 1.16639(1) x 10" 5 /GeV 2 \ 

In order to clarify this strange situation, we make the following observations: 

i) As explained above, both the QED correction C and the electroweak amplitude g 2 / (1— Ar) 
are separately /^-independent quantities. One can certainly evaluate these corrections us- 
ing running couplings. In this case, if the expressions are carried out to all orders, the 
/^-dependence of the couplings is cancelled by the /i-dependence of the radiative correc- 
tions. In actual calculations, due to the necessary truncation of the perturbative series, 
a /i-dependence emerges. However, there is no reason to require that the same scale be 
employed in g 2 /(l — Ar) and in C(/x). Indeed, the natural scale in the former is fi ~ Mz 
and, in the latter, it is somewhere in the range m e < \x < m^, with the optimization meth- 
ods discussed in Section 2 favoring the lower values of \x. For instance, part of C are LB. 
contributions. The photons here carry k 2 = and it is clear that their natural coupling is 
a(m e ) = a, in sharp contrast with a e (Mz). 

ii) Although the choice \i = Mz for the QED correction is not natural, its implementation 
should be done by setting /i = Mz in Eq.([5|), rather than replacing a(m M ) — > a e (Mz) in 
Eq.(p. Using Eq.fl5|) and a e (Mz) to evaluate C(Mz), we find that Gf would be decreased 
by only 0.6 ppm rather than increased by 20 ppm. In particular, we note that changing the 
scale of a(/i) in Eq. ([]]), without modifying the second order coefficient according to Eq.(|), 
is inconsistent at the two-loop level. 
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In summary, our conclusion is that, subject to the neglect of terms of 0(am 2 l /M^ r ), one 
should apply the same QED corrections in the SM as in the Fermi theory. In particular, the 
value Gp = 1.16637(1) x 10~ 5 /GeV 2 , obtained in the Fermi theory, is the one that should 
be applied in the analysis of the SM. 

IV. CONCLUSIONS 

In Section 2, we have applied the traditional FAC, PMS and BLM optimization methods 
to the radiative corrections to the muon lifetime in the Fermi V-A theory. As the corrections 
are convergent, one expects on general grounds the natural mass scale to be in the range 
m e < n < m^. The analysis shows that the FAC and PMS approaches select a mass scale 
very close to m e , while the BLM scheme leads to one very near the geometric average ^/m e m^ 
of the two masses. The FAC and PMS methods provide nearly identical estimations of the 
third order coefficient in the a(m^) expansion, while the BLM estimation is of the same 
sign and differs only by 20% in magnitude. The three optimized expansions give very close 
results, with a maximum absolute difference of 5.4 x 10~ 8 . We use the maximal difference 
between the optimized and a(m M ) expansion, which amounts to be 2.6 x 10 -7 , as an estimate 
of the theoretical error due to the truncation of the perturbative series. This translates into 
an error of 1.3 x 10~ 7 in the determination of Gp, a result very close to estimates obtained 
in Ref. |§ by different considerations. We also find that the expansion in powers of the 
conventional fine structure constant a, which in this case is essentially equivalent to the 
FAC expansion, contains a very small second-order coefficient. In fact, the second-order 
term is about 2900 times smaller than the first. Thus, it turns out that, by a curious 
cancellation of two-loop effects, the original one-loop calculation expressed in terms of 
a, had an error of only 1.4 x 10~ 6 ! Of course, this amusing historical fact could not be 
known before the work of Refs. [[?]|| was carried out. 

In Section 3, we return to the radiative corrections to muon decay in the SM. Using argu- 
ments based on the effective field theory approach, as well as considerations involving higher 
order Feynman diagrams, we show that, subject to the neglect of terms of 0(am^/M^), 
the overall answer factorizes into two separately scale-independent corrections: those of the 
Fermi V-A theory and the electroweak amplitude g 2 /(l — Ar). This important property 
applies to the corrections to the various observables, such as the electron spectrum, electron 
asymmetry in the case of polarized muons, and the muon lifetime. Therefore, if running cou- 
plings are employed, the scales may be chosen judiciously and independently in both factors, 
with fi « Mz being the natural scale in g 2 /(l — Ar), and m e < fi < m M being the logical 
range in the corrections of the Fermi V-A theory. We reach the conclusion that, subject to 
the neglect of terms of C(am^/M^), one should apply the same QED corrections to muon 
decay in the SM as in the Fermi theory. In particular, at variance with a claim presented 
in Ref. ||, we find that the value G F = 1.16637(1) x 10~ 5 /GeV 2 , currently obtained in the 
Fermi theory, is the one that should be applied in the analysis of the SM. 

We conclude our discussion with the following comments: 

i) Although the shifts discussed in Refs. and this paper are numerically very small, 
it is clearly desiderable to evaluate fundamental parameters such as Gf, as accurately as 
possible. 

ii) At present, the Fermi V-A theory is generally viewed, not as independent, but rather 
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as an effective low-energy theory derived from the SM. The arguments presented in this 
paper make clear that the same constant Gp that is precisely determined in the low-energy 
theory is, to a high degree of accuracy, the relevant paramenter in the analysis of the more 
fundamental, underlying gauge theory. 
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